Astronomy & Astrophysics manuscript no. Gj041 


February 2, 2008 


(DOI: will be inserted by hand later) 





The distributed burning regime in Type la supernova models 



F. K. Ropke and W. Hillebrandt 



^1" 

o 
o 

(N 
> 

o 

(N 



> 



o ■ 
Oh; 

6 ' 
a: 



13 



Max-Planck- Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, D-85741 Garching, Germany 
e-mail: fritz ; wfh@mpa-garching.mpg . de 

the date of receipt and acceptance should be inserted later 

Abstract. The deflagration mode of flame propagation in Type la supernova (SN la) models requires a correct description of 
the interaction of the flame with turbulent motions. It is well-known that turbulent combustion proceeds in different regimes. 
For most parts of the deflagration in SNe la the flamelet regime applies. This has been modeled in previous multi-dimensional 
simulations. However, at late stages of the explosion, the flame will propagate in the so-called distributed regime. We investigate 
the effects of this regime on SN la models in a first and simplified approach and show that the trend of effects seems capable of 
curing some problems of current pure deflagration models. 
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1. Introduction 



|Pamkohlerl (fl940) was the first to distinguish between differ- 
ent regimes of turbulent combustion. Since his pioneering work 
progress in the theo retical u nderstanding of the turbulent com- 
bustion process (e.g. Peters, 2000 ) as well as flame ex periments 
in the laboratory (see Abdel-Gaved & Brad levlll98ll for a col- 
lection of data) have confirmed his ideas. He described a regime 
of "large-scale turbulence" where the flame is deformed by in- 
teraction with turbulent motions and a "small-scale turbulence" 
regime where turbulent eddies actually penetrate the internal 
flame structure. In modern perception these regimes are identi- 
fied with the flamelet and th e distributed reg imes, respectively. 

As pointed out by iNiemever & Wooslevl II1997I) . both 
regimes should be reached subsequently in the deflagration 
model of SN la explosions (for a review of models see 
Hillebrandt & Niemeverl |2000j). In a statistical approach the 
distri buted burning regime w as addressed in the context of SNe 
la by iLisewski et alJ (EoOO). Direct numerical simulations of 
flames in degenerate matter h ave recently been a ble to reach the 
distributed burning regime dBell et all |2004b). Nonetheless, 
the distributed burning regime has never been implemented in 
global deflagra t ion models of thermonuclea r supernovae (e.g. 
iReinecke et all l2002l iGamezo et ail l2003h . This was moti- 
vated by the fact that due to the low fuel densities to which the 
regime applies, no additional iron group nuclei are synthesized. 
Moreover, it was assumed that the flame propagation here 
would be too slow to significantly contributing to the energy 
generation. This may be a reasonable approach as long as the 
models aim at first order effects of the explosion characteris- 
tics. However, constant development in modeling and increas- 
ing computational resources facilitate systematic tests of initial 
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parameters (see e.g. Ropke & Hillebrandt! Eo0 4b a) and first 
synthe tic light curves were derived bv ISorokina & Blinnikovl 
(2003). Even synthetic spectra may be calculated from multi- 
dimensional models soon. 

Here the question arises whether the pure deflagration mod- 
els (in which the flame propagates subsonically) in their current 
form are consistent with observational data. The physics input 
of these models may be incomplete. Possible problems arise 
from their low explosion energies which may only reproduce 
the low side of the average observed velocities of the ejecta 
in SNe la. Furthermore, they seem to underproduce interme- 
diate mass elements and leave unburnt material at low veloc- 
ities. Based on their success in one-dimensional parametrized 
models, delayed detonations have b een put forward as a favor- 
able scenario (Gamezo et al., 2004J). In a detonation the flame 
is mediated by shock waves and travels with sound speed. The 
problem with the delayed detonation scenario is that no phys- 
ical mechanism could be identified yet that would trigger the 
transition from a de flagration to a detonation under conditio ns 
of SN la explosions llNiemeverlfl99^lRopke et alll2004albl) . 

With reconsidering the distributed burning regime in the de- 
flagration model of thermonuclear supernovae we propose an 
alternative mechanism that may help to overcome some of the 
difficulties of that model without artificially evoking a detona- 
tion. The advantage of this extended deflagration model is that 
it rests on a sound basis of known physics. If no other effects 
terminate the deflagration phase, burning inevitably enters the 
distributed regime. 

Here, we study the effects of the distributed burning in 
very simplified two-dimensional models. As will be discussed 
below the employed description of this regime is no more 
than a coarse first-order estimate and can only point out the 
trends how it would change the model. Substantially more ef- 
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fort will be needed to implement details of the distributed burn- 
ing regime in global SN la explosion simulations. 

2. Modeling turbulent combustion in SNe la 

The propagation speed of a laminar deflagration flame is de- 
termined by a balance between energy production in the re- 
actio ns and diffu s ive en ergy transport. Based on this simple 
idea, iMikheFsonl (l889) derived an expression for the lami- 
nar burning velocity. This expre s sion w as later corrected by 
IZel'dovich & Frank-Kamenetskvl dl938h taking into account 
the activation energy necessary to induce the reaction. This re- 
sults in a division of the flame structure into a preheat zone and 
an reaction zone. In both theories, the laminar flame speed s\ 
turns out to be proportional to the square root of the diffusivity 
D. 

The wrinkling of the thermonuclear flame in SN la explo- 
sions begins at large scales. An inverse stratification of light 
ashes below dense fuel in the gravitational field of the explod- 
ing white dwarf (WD) star makes burning from its center out- 
ward intrinsically unstable. Buoyancy (Rayleigh-Taylor) insta- 
bilities lead to the formation of burning bubbles that rise into 
the cold fuel. This effect is directly resolved in multidimen- 
sional SN la simulations. Secondary shear (Kelvin-Helmholtz) 
instabilities at the interfaces of the burning bubbles give rise 
to the generation of turbulent motions. The formed turbulent 
eddies decay to smaller scales in a turbulent cascade. 

The eddies of the turbulent cascade interact with the flame 
further wrinkling and stretching it on smaller scales. This leads 
to a rapid increase in the flame surface and results in an en- 
hanced burni ng rate . This effect acts down to the so-called 
Gibson scale (Peters, 2000). Below that scale turbulent veloc- 
ities are so small compared with the laminar burning speed of 
the flame that the flame burns through turbulent eddies before 
they can significantly alter its shape. If the Gibson length is 
large compared with the width of the flame, turbulent burning 
proceeds in the flamelet regime. Here the flame is deformed by 
turbulence, but its internal structure is not affected. 

In the flamelet regime, the effect of surface enlargement on 
numerically unresolved small scales can be compensated by 
attributing an adequate effective turbulent burning velocity s t 
to a (resolved) smoothed flame front. This effective turbulent 
burning velocity has to be determined in a way that the mass 
flux through the smoothed flame equals that through the un- 
smoo thed flame propagating with its laminar speed. lDamko hlei 
dl940t) found that in the flamelet regime the turbulent flame 
velocity completely decouples from the laminar speed and is 
proportional to the turbulent velocity fluctuations v'\ 

s t oc v ' (1) 

Due to the expansion in the explosion process, the 
fuel density drops and th is broadens the flame width (cf. 
iTimmes & Wooslevl Il992l) . At some stage the turbulent mo- 
tions become capable of penetrating the internal flame struc- 
ture. The flame enters the distributed burning regime (see 
lNiemever&KersteinUl997t for an analysis of this transition 
in SNe la). Damktihler's idea to describe this regime was to 
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Fig. 1. Total energies released by the model ignoring dis- 
tributed burning (dashed) and the model including the dis- 
tributed burning phase (solid) . 

assume that the turbulent eddies entering the preheat zone al- 
ter the diffusivity by mixing its thermal structure. Therefore, 
in analogy to the expression for the laminar burning velocity, 
the turbulent burning velocity in the distributed regime should 
scale like 

st/i, ~ (DJDf 12 , 

with D t denoting a turbulent diffusivity. The laminar diffusivity 
D is proportional to product of the flame thickness If and the 
laminar flame speed si. Analogously, D t is set proportional to 
the turbulent velocity fluctuations v' at a certain length scale 
I. For the turbulent flame speed in the distributed regime this 
yields a square-root dependency on v', 

s t ~( Sl v'l/l { ) 1/2 , (2) 

which is confirmed by experimental data. In this regime s t 
scales with the laminar flame speed. Turbulent mixing of the 
internal structure of the flame accelerates its propagation ve- 
locity. 

To model both turbulent burning regimes in numerical 
simulations of SNe l a we appl y the s cheme proposed by 
iReinecke et all dl999l l2002h and iRopkd J2004 . The hydro- 
dynamics is solved for by the piecewise parabolic method 
on a computational grid co-expanding with the WD star and 
the flame is modeled using the level-set approach. In the lat- 
ter, the effective burning velocity of the flame has to be pro- 
vided. To this end v' is derived from a sub-grid scale model 
dNiemever & Hillebrandtlll995h . The maximum value of si and 
expression Q then gives the flame propagation velocity in the 
first stages of the explosion process. 

Previous models stopped burning in the flame once the fuel 
density fell below 10 7 gcirT 3 . To take into account the dis- 
tributed regime we extend this burning and and derive s t form 
Eq. at densities below 10 7 g crrT 3 down to 5 x 10 5 g cm 4 . 
The necessary values for the laminar flame speeds s\ and the 
flame thickness If at low fuel dens ities are provided by fits to 
the data given bv lBell et al.1 d2004albl) . The turbulent velocity 
fluctuations v' on the computational grid scale I are derived 
from the sub-grid scale model. 
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Fig. 2. Snapshots of the models with [1024] 2 cells at t = 10 s. 

The feasibility of modeling the distributed burning regime 
via the level set approach - though in a more elaborate way of 
formulating the flame spee d - has been shown and valid ated 
agains t experiments by iDuchamp de Langeneste & Pitschl 
(2002). 

3. Results of numerical models 

We compare two very simple two-dimensional SN la models 
calculated on a 1024 x 1024 cells grid: one ignoring the dis- 
tributed burning (referred to as "old" model) and one taking 
into account this regime ("new" model ). In both cases the flame 
was ignited in the c3 shape of R einecke et alJ {2002) - a cen- 
tral ignition superposed by sinusoidal perturbations. The initial 
resolution of both models was 1 .975 km. Due to the different 
energy releases the final resolutions were 96 km for the old and 
100 km for the new model. 

The total energy release of both models is plotted in 
Fig. n which also includes results of a resolution study of 
the new model showing convergence for more than [512] 2 



1.00 



x- 0.10 





^c, 16 


"NiV 

i a 

I / 1 1 I 1 1 


// ~ 

/"Mg" 

/ new model, [1024]" cells 

\j old model, [1024] 2 cells 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

t[s] 

Fig. 3. Temporal evolution of the chemical composition in the 
models. 
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Fig. 4. Distribution of unburnt material in velocity space after 
10 s. 



cells. Contrary to more sophisticated three-dimensional models 
the simulations discussed here are based on very simple two- 
dimensional setups and they cannot be expected to reproduce 
the explosion strength of observed SNe la. Nevertheless, trends 
can be inferred even from these simple models. Clearly, the dis- 
tributed regime contributes substantially to the energy genera- 
tion. About 1.1s after ignition the energy generation in the old 
model ceases because the fuel density reaches low values. In 
the new model, however, the burning proceeds for another 0.3 s, 
leading to a ~ 17% increase in the nuclear energy release. In this 
time range, the distributed burning produces large amounts of 
intermediate mass elements. 

This is evident by comparing the snapshots of both simula- 
tions at t — 10 s shown in Fig. [2] They display the position of 
the zero level set as white curve (indicating the position of the 
flame front while burning is active and still providing an ap- 
proximate sep aration betw een fuel and ashes later in the sim- 
ulations; see Ropke, 2004). The distribution of the species is 
color-coded. Our simplified description of the nuclear burning 
(for details see Reinec ke et all |2002) pools the intermediate 
mass nuclei in the representative element "Mg" and represents 
the iron group nuclei by "Ni". The evolution of the chemical 
composition in both models is compared in Fig. [3] 
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Implementing the distributed burning regime into our simu- 
lations we note various important aspects of the results. Due to 
the larger energy release the ejecta reach higher velocities and 
expand faster. This is obvious from Fig. [2 The transition from 
burning in the flamelet regime to the distributed regime takes 
place at low fuel densities (and is implemented in our models 
this way). Therefore the new and old models differ only in the 
production of intermediate mass nuclei. The amount of synthe- 
sized iron group elements is unaffected (cf. Fig.[3J. Moreover, 
Fig. |3 shows that the distribution of "Ni"-rich regions in the 
new model resembles that in the old one. However, in the new 
model these regions are embedded in a layer of intermediate 
mass elements, which also shows a turbulent structure. 

Taking into account the distributed regime, the flame nat- 
urally consumes more fuel. Therefore, the unburnt material in 
the central parts is strongly depleted. This is corroborated by 
Fig. 0] showing an angular average over the distribution of un- 
burnt material in velocity space. No unburnt material is left at 
low velocities. 



4. Conclusions 

By reconsidering the distributed burning regime in global SN la 
models we have proposed a way to alleviate some of the prob- 
lems that might arise when comparing the deflagration model 
with observations. By extending the burning of material to in- 
termediate mass elements it acts against a possible underpro- 
duction of those in previous deflagration models. The addi- 
tional energy released here helps to accelerate the ejecta to 
higher velocities. Near the center the burning consumes now 
most of the fuel leaving no low-velocity carbon and oxygen 
behind. 

The newly implemented burning regime does, however, not 
affect the production of iron group elements since it applies 
to low fuel densities. It was claimed that to reproduce the ob- 
served s pectra a layered com position of the ejecta would be 
needed ( Gam ezo et all Eo04l) . This is, however, not yet con- 
firmed since synthetic spectra derived with multi-dimensional 
radiation transport sche mes are st ill lacking. Recent spectrap- 
olarimetry data (Wang et alll2004l) seem to indicate that a tur- 
bulent structure of the intermediate mass elements is required. 
In any case, the deflagration phase does not remove the mixed 
distribution of species in the ejecta that was present already in 
the old models. 

The results presented here have to be considered prelimi- 
nary for several reasons. Our simulations are performed only in 
two spatial dimensions. The turbulent structure develops differ- 
ently in three-dimensional models resultin g in an overall mor e 
vigorous explosion (e.g. lReinecke et al l l2002llR6pkell2004 . 
An even more important issue is that the modeling of the dis- 
tributed burning regime applied here is very crude and no more 
than a first-order approach. The value of transition density is 
an ad hoc choice. Moreover, Eq. (|2jl is in a strict sense a lim- 
iting case for If » I, which is certainly not yet met at the pre- 
sumed transition density. Therefore our description by simply 
switching from the flamelet regime to the distributed regime 
will overestimate the effects. 



However, the presented study points to the trends that can 
be expected. The distributed burning regime has to be imple- 
mented into global SN la simulations to complete the defla- 
gration model. Important consequences for the energetics and 
the derived synthetic spectra may result from this extension, al- 
though the effects may not be as large as in the crude first mod- 
els presented here. Whether they suffice to remove the current 
problems of the deflagration scenario will be investigated in 
forthcoming studies. For instance, a more general description 
of the turbulent flame propagation velocity ( featuring a sm ooth 
transition between both burning regimes; see Peters, 2000) will 
be tested and also applied to three-dimensional models. 
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